The Origin of Tilted Phase Generation in Systems of Ellipsoidal Molecules with 

Dipolar Interactions 
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We report Monte-Carlo simulation studies of some systems consisting of polar rod-like molecules 
interacting via a pair potential that exhibit liquid crystal phases, attributed with tilt angles of large 
magnitude. For theoretical understanding of the microscopic origin of the tilted phases, different 
systems consisting of prolate ellipsoidal molecules of three different lengths, embedded with two 
symmetrically placed anti-parallel terminal dipoles are considered. We find that the presence of a 
stable tilted phase crucially depends on the molecular elongation which effectively makes dipolar 
separation longer. We observe that in case of mesogens with transverse dipoles the tilt in the layered 
smectic phase gradually increases from zero to a large magnitude as we increase the molecular length. 
However tilt remains weak with molecular elongation for systems with longitudinal dipoles which 
shows a small tilt at shorter lengths. This is the first work determining the combined contribution of 
dipolar separation and orientations in generating biaxial liquid crystal phases with large tilt angles. 
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Polar liquid crystals exhibit a rich variety of tilted lay- 
ered mesophases arising out from the diversity in both 
the inlayer positional arrangements of the molecules and 
in the average molecular orientations with respect to 
the layer normal p], Q • The most common example of 
such a tilted Smectic phase is the Smectic C ( SmC ) 
phase where the center of mass of the tilted molecules 
are randomly distributed in two dimensional fluid like 
layers alongwith a long range orientational ordering of 
molecules over the bulk system. Tilted phases with the 
same point symmetry as SmC phase also exist in other 
orientationally ordered soft matter systems like meso- 
genic polymers, lamellar hp phase and in SmC elastomers 
[3|-[5j. The tilted phases are important not only from 
the fundamental scientific viewpoint but recent discov- 
eries of a large variety of ferroelectric, anti- ferroelectric 
and ferrielectric phase behavior in the tilted phases have 
made them a topic of huge interest from technological 
point of view. These phases are used for developing new 
generation opto-electric devices and various non-display 
applications |6j. 

The origin of tilt in liquid crystalline phases has been 
a topic of much discussion. A large number of theoret- 
ical and experimental studies have been done on tilted 
phases. However the topic remains unsolved in cer- 
tain aspects till today. It is essential to understand the 
structure-property relationship to find the basic interac- 
tions giving rise to tilt in Sm C phase. The microscopic 
origin of the tilted smectic phases are much more com- 
plex than that of the orthogonal smectic phases since the 
existence of a tilted phase is not a favored one because 
packing of tilted rod like molecules in a layer plane needs 
more area than the untilted molecules. This phase can 
be achieved only if there exists additional specific inter- 
actions giving the requisite tilt. 
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In the experimental studies, occurrence of the SmC 
phases has been mostly found in presence of a lateral 
component ofpermanent dipole moment in the organic 
molecules p], [J . As the number of such dipoles increases 
on a molecule, the tendency to form SmC phase also in- 
creases Motivated by the experimental observations, 
Mcmillan gave the mean field theory on the formation of 
SmC phase based on the presence of at least two outboard 
terminal dipoles [8[ . However a freezing of rotation of the 
molecules was resulted during such a SmA to SmC phase 
transition contradicting NMR [9] and neutron scatter- 
ing experiments [lOj which have shown the free rotation 
of molecules in a tilted smectic phase. Wulf attempted 
to relate the formation of SmC phase to the packing re- 
quirements of the zig-zag shaped molecules [11]. However 
lowering of free energy was associated to a freezing of ro- 
tations. A zig-zag model made of seven Lennard Jones 
spheres with two terminal ones at an angle 45° from the 
five in line core, showed a SmC behavior but the tilt ori- 
entation was random and the equilibration process was 
elaborate [12]. A model based on three rigidly linked 
hard spherocylinders arranged in a zig-zag fashion have 
shown the presence of SmC phase [13]. However the sim- 
ple model particles built by assembling in a zig-zag way 
ellipsoidal Gay-Berne particles was unable to show tilted 
phases 1141 . In computer simulation studies by Zannoni 
et.al. [15] , the axial component rather than the trans- 
verse component of lateral dipole moments were found 
to be generating a tilted phase for polar GB molecules. 
In their system, the biaxial order parameters were very 
small and weak non zero tilt was found for the system 
of molecules with lateral dipole moments making an an- 
gle (j) = 0° or 60° with the long axis. Similar results 
were found by Saha et.al. [l6| for a large dipolar oriention 
(j) = 120°. However absence of tilt in these conventional 
model of polar molecules with transverse dipoles, contra- 
dicting experimental results, remains an open problem. 

We explored systems of polar ellipsoidal molecules 
where each molecule is embedded with two terminal anti 
parallel permanent dipole moments placed at equal dis- 
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tance from the center of the molecules to get considerable 
tilt in the smectic phases. In our model, the apolar part 
of the interaction is represented by the Gay-Berne(GB) 
potential which is a modified form of the Lennard- Jones 
potential, considering the anisotropy in interaction. The 
interaction between two such apolar GB ellipsoids i and 
j is given by 

U$ B (Tij,Ui,Uj) = 4e(%,^,%)(pr. 12 - p~ 6 ) 

[Vij - a(fij,Ui,Uj) + a ] . 

where pij = — . Here r^- is the 

separation vector between the center of mass of the 
molecules. The unit vectors U{ and Uj represent the ori- 
entations of the molecules. <to is the minimum separa- 
tion for a side-by-side pair of molecules determining the 
breadth of the molecules. The minimum separation for 
an end-to-end pair of molecules <j e is a measure of the 
length of the molecules. The anisotropic contact dis- 
tance o(fij : Ui : Uj) and the depth of the interaction well 
e(rij,Ui, Uj) depend on the the shape anisotropy parame- 
ter n = (Je/o'o and the energy depth anisotropy parameter 
n = e e /e s which are defined as the ratios of the contact 
distances and energy well depths in the end-to-end and 
side-by-side configurations. The anisotropic contact dis- 
tance varies with n as cr(r^-, Uj) is given by a = 

°L 2 V l + l-x(ui.Uj) J. 

where x — ( ft 2 — 1 ) / ( ft 2 + 1 ) • The anisotropy of the well 
depth e is also controlled by two additional parameters 
ji and v. An explicit description of the GB interaction 
can be found in the original paper [17j. The well depth 
in the cross configuration is written as cq. ctq and eo de- 
fine the length and energy scales. We have used reduced 
units in our calculations by expressing lengths and inter- 
action energies in units of cjq and eo respectively. In the 
present work, we put two point dipole monents on each 
GB molecule at a reduced distance d* = (k — l)/2 along 
the symmetry axis from the center of mass of the uniaxial 
molecules with shape anisotropy We have studied the 
bulk phase behavior for three different values of ft (3, 4 
and 5) keeping the other parameters fixed to their origi- 
nal values n = 5,/i = 2,z/=lin order to investigate the 
effects of varying the dipolar separation to large values. 

The electrostatic interaction energy between two such 
dipolar ellipsoids is given by 

2 *2 
a, (3=1 T a(3 

where r*^(= r*^ — r* a ) are the vectors joining the two 
point dipoles /j,* a and on the molecules i and j at the 
positions r* a = r* ± d*Ui and r*^ = r* ± d*Uj. The re- 
duced dipole moment /i* = y/ \i 2 je^cf^ is chosen /i* =1.0. 
The long range nature of the dipolar interaction is taken 
into account with the reaction field method [18[. The 
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FIG. 1. (color online). Snapshots of the final configurations 
from MC simulations of a system of 1372 GB molecules with 
two terminal longitudinal dipoles for various elongations k : 
(a) Smectic at (k = 3,T* = 1.00, P* = 2.75) with (0) = 
1.6°, (b) Smectic at (k = 4, T* = 1.25, P* = 1.50) with (0) = 
0.6°, (c) Smectic at (« = 5,T* = 1.5, P* = 1.10) with (0) = 
1.6°. 

dipolar part of the total interaction energy including long 
range correction can then be written as 

rrd* = rjd'* _ 2 ( e RF ~ 1) KaV*j(3 

a, 6=1 c 

Where R* = R c /<Jo is the reduced radius of the RF cut- 
off sphere and crf = oo is the dielectric constant of the 
medium. Then the total interaction between two dipolar 
molecules is given by U\f' al — U^ B + Ufj . 
We have performed Monte Carlo (MC) simulations in the 
NPT(isothermal-isobaric) ensemble with periodic bound- 
ary conditions imposed on a system of N=1372 dipo- 
lar molecules. The simulation cell is an orthogonal box 
of dimensions L x1 L yi L z . The dimensions are varied in- 
dependently during simulation so that the system may 
fit itself to its most suitable configuration at each state 
point(P* = Pag/e ,T* = K B T/e ). All the systems 
are prepared initially in a completely disordered isotropic 
phase in a cubic box by melting a crystal structure at 
sufficiently low pressure. We then increase the pressure 
successively by steps of AP* =0.10 or less (near a tran- 
sition). In each case, at a given pressure, the final equili- 
brated configuration obtained from previous lower pres- 
sure is used as the initial configuration. At each state 
point, the system is equilibrated for 3 x 10 5 MC cy- 
cles and 3 x 10 6 MC cycles are used for equilibration 
near a transition. During each MC cycle each molecule 
is randomly displaced and reoriented using metropolis 
criteria where the reorientation moves were performed 
using Barker- Watts technique [18[- One of the three box- 
sides was attempted to change during each MC cycle. 
The acceptance rates of the roto-translational moves of 
molecules and volume moves were adjusted to 40%. 

In order to fully characterize different phases of the 
system various order parameters were computed. The 
average orientational ordering is determined from the 
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FIG. 2. (color online). Snapshots of the final configurations generated by MC simulations of a system of 1372 GB molecules with 
two terminal transverse dipoles for various elongations k : (a) Smectic at (k = 3,T* = 1.00, P* = 1.35) with (0) = 0.25°,(b) 
Tilted Smectic at (k = 4, T* = 1.25, P* = 1.05) with (0) = 6.6°, (c) Tilted Smectic at (k = 5, T* = 1.5, P* = 1.5) with 

(0) = 17°, (d) molecular center of mass positions on a single 2D layer of the tilted phases showing local hexagonal orderings 

(1) = 4, T* = 1.25, P* = 1.05), (2)0 = 5, T* = 1.5, P* = 1.5). 



second-rank tensorial order parameter Q a p defined as 

Q<*p = j? J2iLi(i u ic* u if3- \ & clp) where a, /3 = x,y,z and 
Ui is the molecular end-to-end unit vector of molecule i. 
The nematic order parameter S is given by the largest 
eigenvalue of the ordering tensor Q a p and the corre- 
sponding eigenvector defines the phase director. The 
value of S is close to zero in the isotropic phase and tends 
to 1 in the highly ordered phase. 

We have investigated the smectic structures over var- 
ious elongations k = 3,4 and 5 respectively at fixed 
temperatures T* = 1.00, 1.25 and 1.50. At such tem- 
peratures we observed a jump in S directly from S ~ 
0.10 to S ~ 0.95 indicating a direct isotropic to tilted 
smectic transition. To measure the average tilt angle (6) 
of the phase director about the layer normal, we have 
used the method described in [16] i.e. in smectic phases 
we first find the sets of particles for which first neigh- 
bor distance is < 1.3cr . Each such set forms a differ- 
ent layer. Then we find the normal to each such layer 
by a least square method. The average normal is ob- 
tained by averaging the components over all the layers 
of a single MC configuration. The angle between the 
phase director and the layer normal determines the tilt. 
Then we measure its average value over a number of con- 
figurations. Since all the tilted phases are expected to 
be biaxial we also measure the biaxial order parameter 
(i?2 2) = + cos2 ft) cos 2a cos 27 — cos /3 sin 2a sin 27) 
as described in [19|, where a, /3, 7 are the Euler angles 
giving the orientation of the molecular body set of axes 
w.r.t. the laboratory set of axes. To understand the 
structure of the phases we also calculate the radial dis- 
tribution function g(r) = 47r * 2p (S(r — , where the 
average is taken over all the molecular pairs. 

We have generated well equilibrated tilted smectic 
phases for three different elongations k — 3, 4, 5 with two 
different dipolar orientations = 0° and 90° in order to 
explore separately the roles of a longitudinal component 
and a transverse component in producing tilted phases. 



From the snapshots and various distributions functions as 
described, we analyzed the phase structures. The struc- 
ture of the phases were changed significantly over differ- 
ent elongations. 

In case of mesogens having longitudinal dipoles, the 
GB interaction plays dominant role in the isotropic 
phases and as we increase the pressure along the 
isotherms, the dipolar energy makes larger jump than 
the GB energy in reaching smectic phases. The energy 
distribution in tilted phases show an interesting behav- 
ior over different elongations. For the shortest molecules ( 
k = 3), the GB energy remains stronger than the dipolar 
energy and as we increase elongation they become com- 
parable at k = 4 and 5. The amount of tilt in shortest 
molecules (0) = 1.6° decreases to (6) = 0.6° in the k = 4 
system. The longest (k = 5) mesogens show (0) = 1.6°. 
The snapshots of the phases are shown in Fig.l and the 
corresponding radial distribution functions are shown in 
Fig|3(a)| We see that the amount of interdigitation in 
the smectic phases decreases with elongation showing an 
interdigitated phase at k = 3. 

We now discuss the tilted structures obtained due to 
the effects of two transverse dipoles. In this case, the 
biaxiality comes exclusively from the presence of the 
dipoles. Biaxial smectic phases are found for all the 
three elongations. For n = 3 we obtained an orthogonal 
biaxial phase and for other higher elongations we found 
tilted biaxial phases where the amount of tilt increases 
with k as shown in the Fig. 2. The contribution to 
the total energy is alaways dominated by the dipolar 
interaction in these phases and the dominance increases 
with the elongation. The average tilt order parameter as 
measured gives (0) « 0° for k = 3 and more interestingly 
(0) « 6.6° and 17° for systems with k = 4 and 5 respec- 
tively. We may infer that longer dipolar separation gives 
larger torque arising out from larger dipole moments, 
which is responsible for giving significant change in the 
magnitude of tilt never found before. The values of 
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FIG. 3. (color online) (a) radial distribution functions for systems of molecules with two terminal longitudinal dipoles, (b) radial 
distribution functions for systems of molecules with two terminal transverse dipoles [ The zero of g(r) on the vertical scales 
have been shifted for clarity ] (c) Schematic evolution of the average tilt angle (9) as a function of elongation k for systems of 
molecules with two terminal dipoles . 



the biaxial order parameter in the smectic phases are 
(i?| 2 )=0.87,0.83 and 0.79 respectively for k = 3,4 and 5 
systems. The radial distribution function shows a 
crystalline nature in all three elongations at higher 
pressures as shown in Fig |3(b)| In all the tilted phases, 
the smectic structure were never interdigitated but show 
some difference in local hexagonal ordering as shown 
in Fig. 2(d). Again the tilt angle increased from (6) w 
11° to 17° in smectic phases as we increase the pressure 
P* = 0.45 to P* = 1.5 in systems of longest (k = 5) 
molecules. We noticed that the attempt to increase 
dipole moment by increasing /i* to a value greater than 
1.1 fails due to enhanced probability of dimer formation. 
We have observed the presence of Nematic Phase in 
these systems at different state points as a part of our 
ongoing work. 

Our NPT simulation studies show that the increasing 
dipolar separation as a result of molecular elongation 



can effectively give rise to large tilt in layered liquid 
crystalline phases generated by GB molecules with 
two terminal transverse anti parallel dipole moments 
whereas it has less significant and reverse role in case of 
molecules with two longitudinal dipoles. Experminental 
evidence for tilted smectic phases were earlier reported 
for compounds with two or more lateral dipole moments 
Q, 0. Our study is successful in gaining insights into 
the molecular origin of tilted phases by showing that 
terminal dipoles having longer separation length coupled 
with suitable orientation can bring large tilt to the liquid 
crystal phases. 
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